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Models with essential non-perturbative QCD dynamics and describing Tevatron 
data on high-px charmonium are extrapolated to give predictions of prompt J/ip 
production at the LHC. Differences of up to an order of magnitude occurs. An 
important point is here the treatment of higher order perturbative QCD effects. 



The interplay of hard and soft QCD effects has been demonstrated in 
Tevatron data EJ for high-px J/i 3 , yj' and T , which are up to factors of 50 
above the pQCD prediction in the Colour Singlet Model (CSM) el, where a 
colour singlet cc pair is produced at the parton level. This deficit can be 
explained by letting a fraction of the more abundant colour octet cc pairs 
be transformed into colour singlet through some soft QCD dynamics. This 
has been described m different models: the Colour Qctet (COM) u, Colour 
EvaporationJCEM) Soft Colour Interaction (SCI) □ and Generalized Area 
Law (GAL) □ models, which can be made to fit these Tevatron data. 

In this work we study the extrapolation of CEM, SCI .and GAL models to 
the LHC energy and examine its theoretical uncertainty Dfl. These models are 
based on a similar phenomenological approach, where soft colour interactions 
can change the colour state of a cc pair from octet to singlet. They employ 
the same hard processes to produce a cc pair regardless of its spin state. Of 
importance for high-px J/ip production are NLO tree diagrams with a third 
hard parton that balances the p± of the cc pair. The most important one 
is the process gg — > ccg with a gluon exchange in the t channel, since the 
basic 2^2 process gg — > gg (followed by a gluon splitting g — ► cc) has a 
much larger 0(ct 2 s ) cross section than the LO gg —> cc and qq —> cc production 
processes. Still higher orders than NLO can be important at the Tevatron and 
LHC energies, since many gluons can be emitted and their virtuality need not 
be very large to allow a split into a cc pair. Higher order processes can 
be approximately described by the parton shower approach available in the 
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Figure 1. Distribution in transverse momentum of prompt J/ip as observed by CDF and 
DO 1 in pp interactions at the Tevatron and obtained in the CEM, SCI and GAL models. 

Pythia! Monte Carlo, where in all basic QCD 2^2 processes the incoming 
and outgoing jpaxtons may branch as described by the DGLAP equations. 

In CEM Ux3 the exchange of soft gluons is assumed to randomise the 
colour state, implying a probability 1/9 that a cc pair is colour singlet and 
produces charmonium if its mass is below the threshold for open charm pro- 
duction, m cS < 2mo- The fraction of a specific charmonium state i, relative to 
all charmonia_is_given by a non-perturbative parameter pi ifij/^i — 0.4 — 0.5). 

In SCI it is assumed that colour-anticolour, corresponding to non- 

perturbative gluons, can be exchanged between partons emerging from a hard 
scattering and hadron remnants, leading to different topologies of the con- 
fining colour string-fields and thereby to different hadronic final states. The 
probability to exchange a soft gluon between parton pairs is given by a phe- 
nomenological parameter R. The mapping of cc pairs below the threshold 
for open charm production is based on spin statistics resulting in a fraction 
of a specific quarkonium state i with total angular momentum Ji and main 
quantum number n, given by fa = ^-f^ , where T — (2Jj + l)/nj. 

In GAL i, a generalisation of the area law suppression e~ bA (with A the 
area swept out by the string in energy-momentum space) gives a dynamic 



probability R — i?o(l 



-bAA 



) for two string pieces to interact depending 



on the area difference AA resulting from the changed string topology. This 
favours shorter strings and thereby quarkonium production. The parameters 
Ro and b are obtained from a fit to both HERA and LEP data. 

The comparison of the CEM, SCI and GAL models with the Tevatron 
data is shown in Fig. |l|. As can be seen, all models give a quite decent 
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Figure 2. Differential cross sections in transverse momentum for J/ip in pp collisions at 
t/s = 14 TeV based on the CEM, SCI and GAL models. In all cases the J/ip is required 
to have p± > 10 GeV and decay into fi + fi~ , required to be within the indicated rapidity 
coverage of the Atlas (a) and LHC-b (b) detectors. For comparison, the COM results from 
8 are included in (a). 



description of the data. Although the shape is not perfect in the tail of the 
distribution, it is quite acceptable given the simplicity of the models. The 
overall normalisation is correctly given by the models. For the CEM this is 
obtained by setting pj/^ — 0.43 and the charm quark mass to 1.5 GeV. The 
SCI and GAL models have not been tuned to these data, but the result is 
also sensitive to the charm quark mass (taken as default m c — l-3ix GeV in 
Pythia 5.7). The parton densities used are the default CTEQ3L Il3 for SCI 
and CTEQ4L B for CEM and GAL. One should note that an arbitrary K 
factor is not needed in any of the models, since higher order pQCD processes 
were included through the parton showers. Based on the ability of these 
models to correctly reproduce the Tevatron data and keeping all parameters 
fixed from this comparison, we can now extrapolate them to the LHC energy. 

Results of these models for LHC energies, i.e. y/s = 14 TeV, are shown in 
Fig. U for the acceptance regions of the future LHC experiments Atlas and 
LHC-b. For Atlas , the p± distribution in Fig. ||a is quite similar for the three 
models, although they differ somewhat in the high-pj_ tail. These predictions 
should, however, not be taken as very precise in view of the simplicity of these 
models that attempt to describe unknown non-perturbative QCD phenomena. 
The overall normalisation, which e.g. is sensitive to the value of the charm 
quark mass, should not be considered to be better than within a factor two. 
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For the LHC-b pseudorapidity cuts, SCI is suppressed relative to CEM and 
GAL, showing that the models have a somewhat different rapidity behaviour. 

We can also see in Fig. ||a that the COM result El is significantly lower 
than the other model results. As we have found recently 13, part of this differ- 
ence is because COM only includes cc pairs from the first g — > cc branching 
(namely a NLO approximation), whereas in our calculations with the CEM, 
SCI and GAL models, cc pairs are produced in any branching in the parton 
shower approximation, accounting for all higher orders. For the Tevatron the 
difference between these two approaches is within the precision of data, but 
for LHC the deviation can become as large as one order of magnitude. 

In conclusion, extrapolating models which can describe Tevatron data to 
the LHC energy, we find significant differences in the predicted prompt J/ijj 
cross sections up to almost an order of magnitude. Part of this difference is 
related to the treatment of higher order contributions, which is important in 
order to obtain a good estimate of the correct prompt J /if) production rate. 
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